Elliptic Flow Analysis at RHIC: Fluctuations vs. Non-Flow Effects 
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The cumulant method is applied to study elliptic flow (112) in Au+Au collisions at i/s — 200AGeV, 
with the UrQMD model. In this approach, the true event plane is known and both the non-flow 
effects and event-by-event spatial (e) and V2 fluctuations exist. Qualitatively, the hierarchy of ?Vs 
from two, four and six-particle cumulants is consistent with the STAR data, however, the magnitude 
of V2 in the UrQMD model is only 60% of the data. We find that the four and six-particle cumulants 
are good measures of the real elliptic flow over a wide range of centralities except for the most central 
and very peripheral events. There the cumulant method is affected by the V2 fluctuations. In mid- 
central collisions, the four and six-particle cumulants are shown to give a good estimation of the true 
differential V2, especially at large transverse momentum, where the two-particle cumulant method 
is heavily affected by the non-flow effects. 
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To create extremely hot and dense matter with par- 
tons as its fundamental components - called the Quark- 
Gluon Plasma (QGP) - is a major goal of current and 
future high energy heavy-ion collisions experiments at 
SPS, RHIC and LHC Q. However, due to the complex 
nature of the relativistic nucleus-nucleus reactions, the 
QGP, if it has been created, escapes direct detection. 
Therefore, in order to distinguish the existence and later 
on to investigate the properties of the new state of mat- 
ter, one must find observables which allow to deduce the 
properties of the intermediate (QGP) state from the fi- 
nal state hadrons. Elliptic flow (1*2), which is the sec- 
ond Fourier harmonic |2| in the transverse distribution 
of the emitted particles, is expected to be sensitive to 
the early pressure gradients and therefore the equation 
of state (EOS) of the formed fire-ball in the heavy-ion 
collisions @, 01 • Recent elliptic flow results on Au+Au 
collisions at = 200 AGeV 0, H 0, indeed indicate 
high pressure gradients in the early stage of the reac- 
tion and might therefore hint towards the existence of an 
intermediate QGP state at this energy. 

In principle, the elliptic flow of hadrons at low trans- 
verse momenta (pt) can be related to the degree of ther- 
malization, the viscosity and the EOS of the produced 
matter On the other hand, the elliptic flow of the 

high px particles is related to jet fragmentation and en- 
ergy loss of the primordially produced hard Anti-quark- 
Quark pair when traveling through the hot QCD medium 
[lj| . At low pt, results from most of the RHIC experi- 
ments 0, 0, |El3 indicate a gradual increase of v 2 with 
the increase of pt- This behavior is approximately con- 
sistent with the prediction of relativistic hydrodynamical 
calculation with a first order phase transition to a QGP 
@. When px > 1.5 GeV/c, the vi begins to saturate 
and eventually decreases 0, [ll| . This is a clear signal for 
the breakdown of hydrodynamics at intermediate pt and 
a transition towards jet physics. As shown in [lOj, the 



i>2 at large pt might be a sensitive probe of the initial 
parton density distribution of the Quark-Gluon matter 
produced. Thus an accurate 1*2 measurement might allow 
deeper insights into the bulk properties of the produced 
matter. 

However, an unambiguous experimental measurement 
of the elliptic flow is not a trivial task due to the unknown 
orientation of the reaction plane. Often, experiments use 
the so called reaction plane method |12( to extract the 
magnitude of the elliptic flow. In this method, the re- 
action plane is fixed according to the flow vector of the 
event, then the estimated v-i with respect to the chosen 
reaction plane is corrected for the event plane resolu- 
tion, which accounts for the error in the deduction of 
the reaction plane. The original reaction plane method 
is consistent with the two-particle correlation method 
0, 0, ^3 1 in which v 2 is related to the two particle 
angular difference 4>\ — fa by V2 = \J (cos2(0i — 4>2))- 
However, these two-particle correlations based methods 
might suffer from effects which are not related to the re- 
action plane, these additional contributions are usually 
called non-flow effects 01 j such as the overall transverse 
momentum conservation, small angle azimuthal correla- 
tions due to final state interactions, resonance decays, 
jet production and quantum correlations due to the 
HBT effect ^(|. In order to eliminate the non-flow con- 
tributions to the measured collective flow in the reaction 
plane method, a rapidity gap between the particles used 
to estimate the reaction plane and the measured parti- 
cles is usually introduced. But whether this improve- 
ment works well is still not clear. Recently, the cumulant 
method was proposed fl7| to diminish the non-flow ef- 
fects. The idea of the cumulant method is to extract 
flow with many-particle cumulants, which are the many- 
particle correlations with subtraction of the contributions 
from the correlations due to the lower-order multiplets. 
It is believed that the pure many-particle non-flow corre- 
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lations have much less contributions to the measured flow 
in the many-particle cumulant method. In other words, 
the many-particle cumulant method should be much less 
sensitive to non-flow effects 0] ■ At RHIC ene rgy , the cu- 
mulant method has been applied by STAR 0,0,^3 an d 
PHENIX in the flow analysis of Au+Au collisions. It 
is found that the integral vi from the two-particle correla- 
tions which is denoted as ^{2} is about 15% larger than 
the values from four and six-particle cumulants^ {4} 
and U2{6})- One might attribute the extra v 2 in t>2{2} to 
the non-flow correlations and conclude that the non-flow 
contribution has been successfully eliminated in the re- 
sults with four or six-particle cumulants However, as 
indicated in 0,0], the V2 from many-particle cumulants 
is also affected by the event- by- event v 2 fluctuations. For 
a rough estimation of the fluctuations' contribution to 
the measured elliptic flow, the reader is referred to [T^ |. 
There, the estimation is based on the assumption that 
the i>2 of an event is proportional to initial eccentricity 
of the nucleons or quarks ■ The authors found that 
the difference between «2{2} and t> 2 {4} can also be ex- 
plained by a definite amount of the fluctuations of v% 
which gives a larger v 2 {2} and a smaller «2{4} than the 
exact V2 for the semi-central collisions. Here it is worth 
stressing that this conclusion can be extended to the cu- 
mulant analysis of any order harmonics v n with n > 1, 
i.e. the v n fluctuations result in |t;„{2}| > |f r i{4}| gener- 
ally. Therefore, both non-flow correlations and t>2 fluctu- 
ations can explain the fact that «2{2} > «2{4} at RHIC. 
If the non-flow effects are dominant, the exact V2 is i>2{4} 
or w 2 {6}; if the fluctuations are dominant, the exact i>2 
is (v 2 {2} + u 2 {4})/2 according to 0. This ambiguity 
poses large systematic error to the V2 measurement es- 
pecially at large p t , where the difference between «2{2} 
and i>2{4} is most prominent. So it is essential to make 
clear which effect is dominant in the difference between 
«2{2} and t>2{4} in order to reduce the large systematic 
uncertainty. 

In this article, we use the UrQMD model (v2.2) 00 
to test the robustness of the cumulant method for the 
elliptic flow analysis. The advantages of using a transport 
approach compared to hydrodynamics are immanent: 

• Firstly, transport models do not make any addi- 
tional assumptions on local/global equilibration of 
the matter created during the collisions, but treat 
the non-equilibrium processes directly. 

• Secondly, the present transport approach includes 
most non-flow correlations such as the overall trans- 
verse momentum conservation, small angle az- 
imuthal correlations due to final state interactions, 
resonance decay and jet prodcution, naturally dur- 
ing the systems evolution. 

• Thirdly, the UrQMD model is an event by event 
model, hence it contains the event by event fluctu- 
ations of the elliptic flow. 

Finally, the reaction plane angle <&r is known in the 



model, which allows the direct calculation of the ex- 
act elliptic flow from its basic definition, that is w 2 = 
(cos 2(0 — 3>r)). Therefore, the UrQMD model, even 
if for the time being still underpredicts the integral V2 
in v / s=200AGeV Au-Au collisions at RHIC, is an ideal 
tool to find out whether the V2 fluctuations and non-flow 
effects have large effects on the experimentally used cu- 
mulant method. 

Before the application of the cumulant method, let us 
begin by examining the magnitude of the fluctuations of 
eccentricity and v^. Fig ^ shows a scatter plot of initial 

spatial eccentricity (e = Ppyq4^t ) of the participants as 
a function of impact parameter based on a subset of the 
available UrQMD minimum-bias events. As one can see, 
the eccentricity fluctuations in the model are in magni- 
tude similar to the eccentricity itself. Note that the mag- 
nitude of the fluctuations is quite similar to the estimates 
calculated with a Monte Carlo Glauber model 0]. Due 

to the large event by event fluctuations, (e 2 ) 1 ^ 2 , (e 4 ) 1 ^ 4 

and (e 6 ) 1 are much larger than (e), especially for the 
most central and peripheral events. The eccentricity fluc- 
tuations are supposed to be the main origin of the V2 
fluctuations. Figgis the scatter plot of the event i> 2 av- 
eraged over all particles with \t]\ < 2.5, as a function of 
impact parameter based on the same subset of UrQMD 
minimum-bias events. Like the eccentricity, the event V2 
fluctuations are also of the same magnitude as the V2 

itself. Therefore, (vf) 1 ^ 2 , v y f) 1/ ' 4 an d { v 2) 1 ^ are a l so 
much larger than (i> 2 ), especially in the most central and 
very peripheral centralities where the (i> 2 ) is very small. 
FigGH shows the correlation of i; 2 and eccentricity in the 
semi-central collisions. The fact that the (i> 2 ) is propor- 
tional to eccentricity, shows that the fluctuations in the 
initial conditions are the major origin of V2 fluctuations 
pcj besides the statistical noise [23 due to the limited 
number of particles used to estimate v%. 

The observation of these large fluctuations puts some 
doubt on the accuracy of the experimental methods for 
the extraction of the elliptic flow parameters. Therefore, 
we will now focus on the cumulant method and compare 
the model results (with fluctuations and non-flow effects) 
obtained by different order cumulant methods with the 
exact V2- For the detailed application of the cumulant 
method, the reader is referred to 0|. In our analysis 
we use unit weights in the evaluation of the generating 
function of the cumulants. The parameter r$ is 1.5 as 
usually used in previous analysis. Actually, the detailed 
investigation indicates that the present results are rather 
insensitive to the ro values as pointed out before in |17| . 
For the integral V2 analysis, we use all particles in the 
pseudorapidity region \r)\ < 2.5, and the number of par- 
ticles from event to event fluctuates in each centrality 
bin. We have tested the cumulant method with fixed 
number of particles in each event for the same central- 
ity and find that the results do not change within the 
present statistical error. The centralities in our analysis 
are selected according to the same geometrical fractions 
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FIG. 1: (Color online) Scatter plot of the initial spatial eccen- 

/ 2 \ _ / 2 \ 

tricities (e = <y^) + (^) ) °^ ^ ne participants at different impact 
parameters from the UrQMD model. The black, green, blue 

,1/2 



and grey lines are the average eccentricity (e), ^e 2 )> , ^e 4 )> 



i/4 



and ^e 6 ^ 1 ^ 6 respectively. 




FIG. 2: (Color online) Scatter plot of the event vi averaged 
over all particles within \r)\ < 2.5 at different impact param- 
eters from the UrQMD model. The black, green, blue and 

,1/2 



grey lines are the average elliptic flow {1)2), (wf) > (^2) 



1/4 



and ^wf) 1 ^ 6 respectively. 



of the total cross section (0-5%,5-10%,10-20%,20-30%,30- 
40%,40-50%,50-60%,60-70%) as used by the STAR exper- 
iment Q , however, we use impact parameter cuts instead 
of multiplicity cuts. More than 1.3 ■ 10 6 minimum bias 
events are used in the integral v% analysis. In order to 
increase the statistics at the most peripheral centrality 
bin(60-70%), additional 7 • 10 5 events are added in this 
centrality bin. 

Fig0]shows the calculated integral V2 results as a func- 
tion of centrality. The elliptic flow parameters extracted 
from two-particle cumulant t>2{2} deviate rather strongly 
from the theoretically expected i>2 as obtained from the 
known reaction plane. In fact, at all centralities, i>2{2} 



0.1 



0.05 



- Au+Au 200A GeV (6.4fm<b<8fm) 
'_ average v 2 



-0.05 



-0.2 



0.2 0.4 
eccentricity 



0.6 



FIG. 3: (Color online) Scatter plot of the event vi averaged 
over all particles within \r\\ < 2.5 at different eccentricities 
from the UrQMD model. The red line is the average elliptic 
flow. 



is larger than the exact V2 by 18%. However, for mid- 
central collisions (a/a to t ~ 10 — 50%), the elliptic flow 
from four particle (t>2{4}) and six particle cumulants 
(^2{6}) show almost no difference and both agree well 
with the exact V2- This fact is quite contrary to the 
prediction in |19| . that is, the exact should be in the 
middle of the «2{2} and «2{4} if the differences between 
the cumulant methods are mainly due to V2 fluctuations. 
Therefore, we conclude that for semi-central to semi- 
peripheral centralities the contribution of the V2 fluctua- 
tions to the cumulant results is almost negligible and the 
difference between i>2{2} and i>2{4} or i>2{6}, is mainly 
due to non-flow effects in the UrQMD model. The sim- 
ilar conclusion is also made in |24j . where the directed 
flow vi in 1.69A GeV Ru+Ru collisions at SIS/GSI is 
measured with the cumulant method. It is found that 
|«i{2}| is less than |^i{4}| which is also contrary to the 
predictions based on flow fluctuations. 

In spite of the good agreement of ^{4} or V2{6} with 
the exact i>2 in semi-central bins, from Fig0J we have 
also seen that both «2{4} and V2{6} do not agree with 
the exact V2 in the most central and the very peripheral 
bins. This means at central and very peripheral colli- 
sions, the i>2 fluctuations indeed play an important role 
as indicated in |l9j. In the peripheral bins the higher 
order cumulants give larger V2 than the exact one. In the 
most central bin, the i>2{4} is smaller and even becomes 
complex (not shown in FigQJ due to the fluctuations, 
while the «2{6} is slightly larger than the exact V2- These 
findings are qualitatively consistent with previous results 
within a simplified Monte-Carlo Glauber treatment [Tflj . 

In order to estimate how sensitive the cumulant 
method is to impact parameter fluctuations in a central- 
ity bin, we also performed the cumulant analysis in en- 
larged centrality bins. The pink (grey) points in Fig0] 
show the results for the enlarged bins(0-10%, 5-20%, 10- 
30%, 20-40%, 30-50%, 40-60%, 50-70%). One can see 
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that the vi values from any order cumulants are still in 
line with the corresponding t> 2 results from the original 
bins although the impact parameter fluctuations in the 
enlarged bins are larger than those in the original (nar- 
rower) centrality bins. Thus, the main contribution to 
the V2 fluctuations in the original centrality bins should 
be due to V2 fluctuations at the same impact parame- 
ter, e.g. due to the spatial eccentricity fluctuations from 
event to event and not due to impact parameter fluctua- 
tions. 

While the total elliptic flow values extracted from the 
calculation are lower than the experimental results, the 
relations between «2{2}, i>2{4} and «2{6} are similar to 
the results reported by the STAR collaboration at RHIC. 
As shown in FigEJA), open symbols denote the calcula- 
tion, while full symbols show the STAR data on the ratios 
v 2 {2}/u2{4} and f2{6}/i>2{4} for comparison. The good 
agreement between UrQMD results and the data may 
indicate that the mechanism which accounts for the dif- 
ferences between t>2{2} and i>2{4} or 1)2(6} is the same. 
In Fig[5{B), we show the §2 factor from the UrQMD 
model. The (72 factor, is defined in experiments as |25| 
92 = N- (w2{2} 2 — 1>2{4} 2 ), where N is the event multiplic- 
ity (for our analysis) or the number of wounded nuclcons 
(for the STAR data) which should be approximately pro- 
portional to the multiplicity. The 52 is a measure of the 
non-flow effects and should be independent of the central- 
ity as originally suggested by [25|. However, the STAR 
Q and SPS [26| data show that with the increase of the 
impact parameter, the gi will decrease by about a fac- 
tor of 3. This decrease of the observed gi is consistent 
with the results based on the eccentricity (or ^2) fluctua- 
tions |l9j, which seems confirm the conjecture in [25j, i.e. 
the V2 fluctuations account for the variance of 52 with 
centrality. As we can see in Fig|5jB), the gi from the 
UrQMD model (the blank square) also has similar shape 
as the data (please note that gi from UrQMD has been 
rescaled by a factor 0.186 to compare to the 200AGeV 
STAR data, since the magnitude of the V2 is too small 
and N in the g2 defintion is also different). However, 
the decrease of <?2 in the UrQMD model is not mainly 
from i>2 fluctuations. Since the exact V2 is known in the 
UrQMD model, the exact two-particle non-flow correla- 
tions (or the exact 52), i-e- A • («2{2} 2 — t>f), is also 
known. In FigJSjB), the solid line is the exact g2 from 
the UrQMD model. Apparently, the exact §2 decreases 
towards peripheral centralities too. Therefore, although 
the V2 fluctuations indeed make the experimentally mea- 
sured <?2 decrease slightly faster, the non-flow correlations 
themselves change with the centralities too. 

Recently, to overcome the experimental limitations in 
the V2 measurement with the reaction plane method, the 
STAR experiment has upgraded its set-up. The Shower 
Max detector of the Zero Degree Calorimctcrs(ZDC- 
SMD) has been added to reconstruct the reaction plane 
with the sideward deflection (bounce-off) of the spectator 
neutrons. The non-flow effects are supposed to be min- 
imal, because the spectator neutrons barely participate 



in the complicated final state rescattering. The STAR 
preliminary results |27| for the measured i>2 with respect 
to this reaction plane is denoted as w 2 {ZDC-SMD}. The 
reported «2{ZDC-SMD} agrees well with i>2{4} in the 
mid-central collisions (10-50%). v 2 {ZDC-SMD} is larger 
than i>2{4} in the most central bins(0-5% and 5-10%) and 
smaller than the «2{4} in the very peripheral bins (larger 
than 50%). The relation between ^{ZDC-SMD} and 
V2{4} is similar with those between exact V2 and «2{4} 
from UrQMD. This similarity, on the one hand, confirms 
that the mechanism which affects the cumulant method is 
indeed the same as that in UrQMD; on the other hand, it 
supports that the flow measurement with the ZDC-SMD 
is not disturbed by non-flow effects or flow fluctuations. 
Therefore we want to advocate the ZDC-SMD method 
for the flow analysis, because it allows to extract very 
reliable flow results over the whole centrality. 
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FIG. 4: (Color online) The integral i>2 results(i>2{2},U2{4} and 
U2{6}) from the cumulant method are compared to the exact 
V2 in different centrality bins. The pink (grey) points are the 
corresponding results from the enlarged centrality bins which 
merge two of the original bins. 

Let us now turn to the study of the the differential 
r>2- In the cumulant method, the differential V2 in one pt 
or rapidity bin is estimated with the cumulants between 
the particles in this bin and those in one common "pool" . 
The average V2 of the particles in the "pool" should be 
known from the integral flow analysis. For the following 
differential V2 analysis, we always use all the particles 
within < 2.5 as the "pool". One should also notice 
that the non-flow correlations which affect the differential 
flow analysis will be that between the particles in the 
chosen bin and those in the "pool" . 

Firstly, let us explore the px dependence of i>2- Here we 
use more than 6 TO 5 semi-central events (with impact pa- 
rameters from 6.7 to 8.3 fm corresponding to about 20% 
to 30% of the total cross section). From the above re- 
sults on the integral V2 , we know that both four and six- 
particle cumulants produce almost the exact V2 in this 
centrality bin, but it is still necessary to see whether 
the different cumulant method produce the differential 



5 



as 



1.6 
1.5 
1.4 
1.3 
1.2 
1.1 t- 
1 I 
0.9 



0.2 



0.1 



(A) 



UrQMD v 2 {2}/v 2 {4} 
UrQMD v 2 {6}/v 2 {4} 
STAR v 2 {2}/v 2 {4} 
STAR v 2 {6}/v 2 {4} 



+-: 



-ir- 



IB) 



STAR 200 
STAR 130 
UrQMD 200 




10 20 30 40 50 
% Most Central 



70 



FIG. 5: (A) The ratios v 2 {2}/v 2 {4} and v 2 J6}/v 2 {4} from 
UrQMD are compared to the STAR data Q{. (B) The g 2 
factors from the UrQMD model are compared to the STAR 
data, the solid line is the exact g 2 , i.e. N- (v 2 {2} 2 — v 2 ) (please 
see the text for details). Note that g 2 's from UrQMD have 
been scaled down by a factor 0.186 



V2 correctly. Especially at large transverse momenta (px) 
non-flow contributions are expected to be large and might 
influence the results obtained by the cumulant method. 
FigHj] shows the calculations for the vi of particles within 
\r\\ < 2.5 as a function of px- At low px, the exact vi 
increases with the increase of pt and reaches a maximum 
at about 2.5 GeV/c, then drops down with a further in- 
crease of px- In contrast, the vi from two-particle cu- 
mulant i>2{2} increases also at low px, but stays roughly 
constant at large px, in addition it is always higher than 
the exact vi- The saturation of f2{2} is consistent with 
STAR's «2{2} results This strong deviations point 
towards substantial contributions from non-flow effects 
in the two-particle cumulant method. The higher or- 
der cumulants do a much better job in reproducing the 
exact V2- Here, the difference between «2{4} and the ex- 
act V2 is much smaller especially at large px- However, 
D2{4} is still larger than the exact v%, indicating that even 
four-particle cumulants are not free from non-flow distur- 
bances. When we go to the six-particle cumulant results 
f2{6}, we get good agreement with the exact V2 in the 
whole px range within the statistical error. This shows 
that the non-flow effects have been completely eliminated 
in v 2 {6}. 

Finally, we will study the pseudo-rapidity (rf) depen- 
dence of i>2 with the cumulant method using the same set 
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FIG. 6: v 2 (pt) in the semi-central collisions: results from the 
cumulant method are compared to the exact v 2 



of semi-central events as for transverse momentum anal- 
ysis. It is usually expected that at large 77, the non-flow 
effects are less important than at mid-rapidity because 
of the larger rapidity gap between the particles in the 
rapidity bin and the "pool" particles. So the difference 
between 7^2} and t>2{4} might be smaller at large 77 
compared to midrapidity. FigUA) shows the results on 
112(77) obtained from the different methods. Indeed one 
observes that at large 77, «2{2}, «2{4} and «2{6} are al- 
most similar and they all agree well with the exact x>2- 
This is in line with the STAR results on the ^2(77) also 
indicating agreement between «2{2} and 7^4} at large 77 
0. However, the smaller difference between the t^'s from 
any-order cumulants at larger rapidity must not be taken 
as a sign that the non-flow effects are less important at 
larger rapidities, because the V2 itself decreases towards 
large rapidity. To demonstrate this, Fig0B) shows the 
ratios of V2 {n} over the exact V2- One observes that the 
ratios are roughly independent of the rapidity. There- 
fore, the non-flow effects at forward rapidity might be as 
important as those at mid-rapidity. 

After the test of the cumulant method, we know that 
V2 from the two-particle cumulant is heavily affected by 
the non-flow effects. While the still remaining question 
is what is the main origin of these non-flow correlations. 
Firstly, in the UrQMD model, there is no correlations 
due to the HBT effect. The correlations due to the over- 
all transverse momentum conservation are also found to 
be less important for the t;2 measurement ^4j . The short 
range direct correlations due to resonance decay, jet pro- 
duction and the last elastic or inelastic interaction before 
freeze-out can be removed by selecting the particles in 
the analysis preferentially. That is to say, only one of 
the final state particles that are directly from the above 
processes is selected in the analysis. The integral V2 {2} 
and ^{4} of the selected particles in the semi-central 
bin (20-30%) are almost the same as those of the original 
analysis. Therefore, the short range direct correlations 
are not the main non-flow correlations at least for this 
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FIG. 7: (Color online) (A) 1)2(77) in semi-central Au+Au colli- 
sions at y/s = 200 AGeV. Results from the cumulant method 
are compared to the exact V2- (B) Ratios of V2{n} over the 
exact V2- 



centrality. It should be stressed that the jet processes 
happen at the very beginning of the event and most res- 
onances decay before the freeze-out of the system. There- 
fore, most of the daughter particles directly from jets or 
resonances will interact with other particles in the sys- 
tem. And the correlations between these daughter par- 
ticles could be changed by the final state interactions. 
As has been shown in j2^, the final state rescatterings, 
or more exactly, the transverse radial expansion can con- 
tribute to the non-flow correlations between these daugh- 
ter particles due to the fact that the particles from the 
same NN collsions (or jets and resonances) are pushed in 
the same direction radially in the transverse plane. Ref. 
pEjj predicts that, this kind of non-flow correlations will 
depend on the centrality (following the development of 
radial flow). It is indeed observed both in the UrQMD 
results and in the data that 52 factor increases towards 
central collisions, which shows the non-flow correlations 
between the daughter particles of jets and resonances due 
to transverse radial expansion could be very important 
or even dominant for semi-central to central collisions. 

In summary, we have applied the cumulant method 
to analyze the 1*2 of the Au+Au reactions at y/s = 
200 AGeV within the UrQMD model. On the integral 
i>2 analysis, we reproduce the hierarchy of i>2{2}, «2{4} 
and «2{6} observed by the STAR experiment even if the 
V2 from UrQMD is only about 60% of the data. From 
the comparisons of the cumulant results to the exact V2, 



we found that v-i fluctuations affect the results from the 
cumulant method in the most central and very peripheral 
collisions. However, this effect is almost negligible over 
a wide range of the mid-central collisions (about 10-50% 
of the total cross section). 

While the two-particle cumulant results are heavily af- 
fected by non-flow effects, non-flow effects can indeed be 
nearly eliminated using four and six-particle cumulants. 
The similarity between STAR data and UrQMD results 
shows that the new flow measurements at STAR (using 
the ZDC-SMD detector) are a good way to obtain V2 val- 
ues which are not disturbed by the non-flow effects and 
V2 fluctuations over the whole centrality range. 

For the differential V2 analysis, the two-particle cumu- 
lant method gives a nearly saturated V2 at large px in 
stark contrast to the exact V2 that drops down rapidly 
at high pt- The V2 from four and six-particle cumulants 
agree well with the exact «2 especially at large px- How- 
ever, there are still some non-flow contributions left in 
the four-particle cumulant method so that the 1^2 {4} is 
always a little (about 4%) larger than the exact V2- Fi- 
nally, we point out that in the present model the non-flow 
effects at forward rapidity might be as important as those 
at mid-rapidity. 

As for the origin of the non-flow correlations, we find 
that the short-range direct correlations due to resonance 
decay, jet production and the last final state interactions 
before freeze-out are not important for the semi-central 
to central Au-Au collsions at RHIC. While the indirect 
correlations between the daughter particles of jets and 
resonances could be the major non-flow effects for these 
collisions, especially these correlations could be substan- 
tionally enhanced by the transverse radial expansion of 
the medium |2S| . 

A final remark. As shown above, the many-particle 
cumulant method V2{n > 4} allows for a good estima- 
tion of the exact V2 in the mid-central collisions, thus 
one may justify other analysis methods by comparing 
their results with the cumulant method results. For in- 
stance, the PHENIX reaction plane method [f| seems 
also to suffer from non-flow effects because it gives the 
same V2(pt) results as the two-particle cumulant method 
which is heavily affected by the non-flow effects as dis- 
cussed above. The STAR ZDC-SMD method seems to 
give good estimates of the integral flow. But further 
comparisons with the cumulant method on the differen- 
tial flow are necessary to fully justify its application in 
the differential flow analysis. 
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